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MOTE  ON  SHEAR  INTERACTION  BETWEEN  TWO  FIBERS 


Introduction 

The  shear  lag  approximation  for  stress  distributions  in  damaged  composites 
employed  by  Hedgepeth  and  Van  Dyke  [1],  and  later  by  others,  employs  an 
effective  shear  stiffness  Gh  to  describe  the  role  of  the  matrix  in  transferring 
direct  stress  from  a  fiber  to  a  neighboring  fiber. 

In  their  treatment  it  was  assumed  that  the  force  per  unit  length 
transferred  between  fibers  1  and  2  Is  given  by 

F  ’  Af  ^2  *  91  ^2  *  *l)/d 

where  Af  is  fiber  area,  a  is  fiber  stress,  G  is  the  matrix  shear  modulus,  d 
Is  distance  between  fiber  centers,  and  w2  -  w1  is  the  relative  axial  displace¬ 
ment.  Shear  lag  was  originally  devised  as  a  technique  for  obtaining  stress 
distributions  In  aircraft  sheet-stiffener  combinations.  In  that  context  d 
is  the  distance  between  rivet  lines  and  h  is  sheet  thickness.  In  extending 
the  analysis  to  a  3D  composite  d  remains  well-defined,  but  loses  physical 
significance,  since  (w2  -  Wj)/d  no  longer  serves  as  a  direct  measure  of  shear 
strain.  But  h  is  undefined,  and  loses  physical  significance  as  well. 

Hedgepeth  and  Van  Dyke  said  nothing  about  how  to  evaluate  h.  For  their 
purposes  its  magnitude  did  not  matter.  They  only  used  shear  lag  theory  to 
evaluate  stress  concentration  factors  In  cracked  uni  directionally  reinforced 
composites  having  fibers  in  square  or  hexagonal  arrays,  and  they  found  that 


(1) 


2 


for  this  problem  the  results  are  Independent  of  the  value  of  h.  However,  the  main 
reason  for  wanting  to  know  stress  concentration  factors  Is  to  evaluate  the 
probability  of  crack  extension.  This  requires  knowing  not  only  the  stress 
concentration  factor  but  also  the  entire  stress  distribution  in  fibers  borderinga 
crack,  and  stress  distribution  does  depend  on  h.  Thus  the  theoretical  analysis  of 
damage  accumulation  requires  being  able  to  evaluate  h.  While  lack  of  know¬ 
ledge  of  h  is  most  serious  in  the  case  of  3D  composites,  it  is  also  a  handi¬ 
cap  in  the  study  of  2D  composites. 

Different  authors  have  used  different  tactics  in  studying  damage  accumulation 
and  comparing  theory  and  experiment  without  knowing  h.  For  instance,  Rosen  and 
Zweben  [2,3]  assumed  the  stress  concentration  factor  is  constant  over  virtually 
the  entire  overload  region  in  a  fiber  neighboring  a  crack.  In  calculating 
the  extent  of  this  region  they  used  an  analysis  in  which  a  broken  fiber  is 
embedded  in  a  cylinder  having  the  properties  of  the  matrix  and  both  are 
surrounded  by  an  infinite  orthotropic  elastic  medium  having  the  average 
properties  of  the  composite  [4],  Goree  and  Gross  found  the  stress  distri¬ 
bution  in  neighboring  fibers  using  a  modified  form  of  shear  lag,  and  in  an 
early  study  [5]  assumed  in  effect  that  h  and  d  are  equal.  This  cannot  be  true, 
since  it  implies  that  the  load  transferred  between  fibers  is  independent  of 
their  separation,  whereas  it  is  easily  shown  that  the  load  transfer  approaches 
infinity  as  the  distance  of  closest  approach  goes  to  zero.  In  later  studies 
Goree  and  colleagues  have  evaluated  the  effective  shear  transfer  stiffness 
|  for  a  given  composite  by  comparing  the  theoretically  calculated  and  the 

j  experimentally  measured  crack  opening  displacement  [6,7].  Batdorf  and 

t 

Ghaffarian  [8]  did  not  use  h  at  all,  but  found  the  effectiv*  size  of  the 
overloaded  region  in  fibers  adjacent  to  cracks  in  a  given  composite  by 
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comparing  4  theory  for  ultimate  strength  with  experimental  data.  The  last  two 
approaches  can  only  be  applied  to  composites  that  have  already  been  fabricated 
and  with  which  the  appropriate  experiments  have  been  carried  out. 

It  seems  likely  that  In  the  30  case,  h  depends  on  fiber  size,  fiber 
separation,  and  composite  geometry  [square,  hexagonal  or  some  other  array]. 

An  approximate  value  for  it  could  probably  be  obtained  using  the  methods  of 
[4].  A  more  accurate  value  can  be  found  experimentally  using  an  electric 
analogue  [9].  A  thorough  evaluation  using  the  analogue  will  require  some 
time,  however.  The  present  paper  gives  an  exact  analytical  treatment  of  h 
fbr  the  limiting  case  of  a  single  pair  of  fibers.  This  will  provide  a  standard 
against  which  to  measure  the  accuracy  of  the  electric  analogue  approach  and  may 
shed  some  light  on  the  value  of  h  for  a  composite. 

Theory 

The  theory  for  shear  transfer  in  the  case  of  two  Infinitely  long  cylinders 
bonded  to  an  infinite  elastic  matrix  closely  parallels  problems  with  a 
similar  geometry  in  such  fields  as  electrostatics,  fluid  mechanics,  potential 
theory,  and  complex  variables.  Those  already  familiar  with  such  problems 
may  prefer  to  go  directly  to  (4).  Here  we  offer  an  elementary  solution 
specifically  relating  to  the  problem  In  hand. 

Consider  a  line  in  an  infinite  elastic  medium.  The  line  passes  through 
the  point  (-a,o)  in  the  x-y  plane  and  extends  to  ±  ®  in  the  z-dlrection. 

If  It  Is  loaded  with  an  axial  force  F'  per  unit  length,  there  will  be  a  dis¬ 
placement  w  (x,y)  In  the  surrounding  elastic  medium.  Since  the  axial  load 
transmitted  from  the  Inside  to  the  outside  of  any  coaxial  cylindrical  shell 
must  be  Independent  of  radius, 

F'  ■  2m  *-2trr  G  fa 

The  negative  sign  Is  used  so  that  w  and  F'  will  have  the  same  sign  convention. 


(2) 
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From  this 

*(r)  ="Hg  £n  ^ 

Here  r  is  the  distance  of  the  field  point  from  the  line  load,  and  G  is  the 
shear  modulus  of  the  medium.  RQ  is  a  large  distance,  approaching  infinity, 
at  which  w  is  essentially  zero. 

If  an  equal  and  opposite  line  load  is  applied  at  (a,o)  the  resulting 
axial  displacement  at  (x,y)  will  be 
w(x,y)  =-^  An  Q  ■  In 

*  **  F'  £n  (r,/rj  (4) 

2tFg 

Here  r^  and  r?  are  the  distances  of  the  field  point  (x,y)  from  the  load 
points  (-a,o)  and  (a,o)  respectively. 

The  loci 

rl^r2  *  c  =  constant  (5) 

are  cylinders,  the  traces  of  which  In  the  x,y  plane  are  as  shown  in  Figure  1. 

On  each  cylindrical  surface  the  displacement  w  is  constant.  Accordingly 
the  elastic  materials  inside  the  cylinder  can  be  replaced  by  a  rigid  rod 
without  altering  the  stress  distribution  outside  of  the  cylinder.  The 
cylinders 


rl/r2  *  C 

(6a) 

1 

and 

r?/r 1  »  c 

(6b) 

are  mirror  Images  of  each  other  with  respect  to  the  xz  plane.  They  will 

K 

sncnSja 
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represent  for  our  purposes  two  equal  and  parallel  rigid  rods,  mutually 
displaced  axially  a  distance  such  that  the  load  transferred  per  unit  length 
is  F',  namely. 


“2  -  »1  *  IS  tn  c 


F‘  *  irG  (w?  -  Wi)/tn  c 


We  now  seek  to  use  the  above  analysis  to  evaluate  h  in  Hedgepeth's 
equation  (1).  We  note  that  in  terms  of  the  set-up  In  Figure  1,  the  radius 
of  a  circle  described  by  (5)  is 


'(ftr) 


Two  such  circles  have  their  centers  at 


c2*l  , 

*  ■  ‘TTT 


Thus  the  distance  between  centers  is 

.  _  o(c2  +  1)  . 


(c*  -  1) 
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Comparing  (1)  and  (7b) 
h/d  «  rr/£n  c 

The  distance  of  the  closest  approach  can  be  shown  to  be  given  by 


(11) 


(12) 


Thus 


»  (g-'-U.  .%li2  03) 

(cZ  +  1)  (c  +  1)  c  +1 

(11  )  and  (13}  may  be  regarded  as  parametric  equations  relating  h/d  to  s/d. 
Alternatively  (13)  may  be  solved  for  c  as  a  function  of  s/d  and  the  result 
substituted  in  (11),  with  the  somewhat  cumbersome  result 

h/d  =  n/ln  [(1  t  /|.(i_s/d)2  )/(l-s/d)]  (14) 


This  relationship  is  shown  graphically  in  Figure  2. 

Discussion 

It  is  obvious  from  Fig.  2  that  h/d  is  not  a  constant,  but  varies 
with  the  ratio  s/d.  As  s  -*•  0,  h/d  »  as  it  should,  and  for  s/d  -  1 , 
h/d  ■+  0.  For  s/d  *  0.2,  the  value  for  a  composite  having  a  fiber 
volume  fraction  of  0.5,  h/d  5  4.5. 

Equation  (14)  represents  an  exact  solution  of  the  2-flber  problem. 
The  approach  employed  cannot  be  extended  to  a  larger  number  of  circular 
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fibers  because  for  that  case  the  w  «  constant  contours  are  no  longer 
circular.  The  electric  analogue  approach  can  be  used  to  obtain  very 
accurate  solutions  for  arbitrarily  chosen  values  of  the  geometric  para¬ 
meters,  but  of  course  does  not  yield  an  analytic  solution  such  as  that 
given  here  for  the  two- fiber  case. 
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